Determination of the chiral condensate from 2+1-fiavor lattice QCD 
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We perform a precise calculation of the chiral condensate in QCD using lattice QCD with 2+1 
flavors of dynamical overlap quarks. Up and down quark masses cover a range between 3 and 
100 MeV on a 16 3 x 48 lattice at a lattice spacing ~ 0.11 fm. At the lightest sea quark mass, 
the finite volume system on the lattice is in the e regime. By matching the low- lying eigenvalue 
spectrum of the Dirac operator with the prediction of chiral perturbation theory at the next-to- 
leading order, we determine the chiral condensate in 2+1-flavor QCD with strange quark mass fixed 



at its physical value as £ MS (2 GeV) 
systematic, respectively. 



[242(04X1}?) 



MeV] where the errors are statistical and 



Spontaneous breaking of chiral symmetry is one of the 
most fundamental properties of quantum chromodynam- 
ics (QCD), as it produces the bulk of the hadron masses. 
The symmetry breaking is indicated by a nonzero value 
of the chiral condensate £, which is an expectation value 
of the scalar density operator qq. Despite its importance, 
calculation of £ remains a significant challenge, even us- 
ing the numerical simulation of QCD on the lattice, due 
to both ultraviolet and infrared problems. 

On the ultraviolet side, an additive renormalization 
of the scalar operator diverges as ~ 1/a 3 as the lattice 
spacing a decreases, when the chiral symmetry is vio- 
lated. Even with exact chiral symmetry, there exists a 
quadratic divergence proportional to the quark mass. On 
the infrared side, since spontaneous symmetry breaking 
does not occur at finite volume, the infinite volume limit 
has to be taken before going to the massless limit. There- 
fore, careful study of the scaling in the chiral and infinite 
volume limits is crucial to determine £. 

Our previous work 0, [2] opened a new possibility to 
overcome these difficulties by performing a lattice QCD 
employing the overlap fermion formulation 0, Q , which 
preserves exact chiral symmetry at finite lattice spacings. 
The ultraviolet problem is avoided by using the spectrum 
of low-lying fermion modes. According to the Banks- 
Casher relation [5], the spectral density p(X) of the Dirac 
operator at A = is related to the chiral condensate as 
£ = np(0). At a large but finite volume V, chiral pertur- 
bation theory (ChPT) can be used to predict the volume 
scaling of the near-zero modes, which is also equivalently 
described by the chiral random matrix theory By 
matching the theoretical prediction with the lattice data, 
the chiral condensate E was determined at the leading 
order (LO) in the e expansion (See also (lOj). 

This letter extends the previous work in several direc- 
tions: (i) Based on a new ChPT calculation by Damgaard 
and Fukaya (llT |. which is valid in the conventional p 
regime as well as in the e regime, we use the lattice data 
at several values of sea quark masses, (ii) The new for- 



mula consistently treats the next-to-leading order (NLO) 
effects in the p expansion and thus, the result of £ has 
the NLO accuracy (A similar NLO analysis of the lattice 
data taken with the Wilson fermion in the p regime has 
been done recently [12j|.). (hi) The lattice data are newly 
generated including the effect of strange quark, so that 
the result corresponds QCD in nature, (iv) The finite 
volume scaling is confirmed using two volumes 16 3 x 48 
and 24 3 x 48. With these new developments, the deter- 
mination of E is made more precise and reliable. 

The spectral density at a given topological charge Q is 
calculated within ChPT at NLO as [TO] 

Pq(X) = £ eff/ Sn(AE eff V, {m sea 'E cS V}) + p p (\, {m sea }), 

' ^ (') 

for an eigenvalue A of the Dirac operator. Assuming the 

analyticity, Pq{X) is obtained through the real part of 
the chiral condensate with a valence quark mass equal to 
an imaginary value iX. Here E e g is an "effective" chiral 
condensate of which definition is given below. 

The spectrum of the near-zero quark modes (A ~ 
1/EV) is mainly affected by the zero-momentum pion 
modes. In fact, the first term in ([T]) is the same as 
the spectral density at the leading order of the e ex- 
pansion expressed as a function of dimensionless 
combinations XE e sV and {m sea 'E c ffV} — {mi'E e sV, • • • , 
mjVfEoffy} expressed by 



Pq(C,{/W}) = C 2 
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with Nf x Nf matrix A and (N f + 2) x (N f + 
2) matrix B defined by Aij — p{ 1 lQ+j-\{p,i) and 

By = e 2 J Q i.i 2(C)- B 2j = C^JQ+.-iiO, Bit = 
{-Pi-2) 3 1 lQ+j-i(Pi-2) (i ^ 1,2), respectively (J fe 's 
and Ii's denote the (modified) Bessel functions.). The 
phase factor C2 is 1 for Nf = 2 or 3. 

The second term in (p} is the NLO correction seen 
in the ordinary p expansion With the meson mass 
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Mf- = (nii + mj^/F 2 , which is made of either sea quark 
(/) or valence quark (v), it is expressed as 



ttF 2 



Re 



^T(A(M 2 v )-A(M 2 ff /2)) 
f 



-(G(Ml)-G(O)) 



(3) 



The function A(M 2 ) contains the chiral logarithm, 
A(M 2 ) = ^ In If- + gi(M 2 ), with gi(M 2 ) represent- 
ing the finite volume effect [13] ■ The subtraction scale 
(i su b is set at 770 MeV in this work. The other function 
G(M 2 ) has a double-pole contribution due to the par- 
tial quenching. The explicit forms of gi(M 2 ) and G(M 2 ) 
are given in [ll|[In this Letter, we use a simplified nota- 
tion. G(M 2 ) corresponds to G(0,M 2 ,M 2 ) in 1 1 jj . . The 
effective chiral condensate £ e ff in (fTJ) is given by 



-(5(0) -leLj^ 




where Lg (renormalized at /^ s «b) is one of the low-energy 
constants at NLO (la ). 

Numerical simulations of lattice QCD are performed 
using the Iwasaki gauge action at j3 = 2.3 including 2+1 
flavors of dynamical overlap quarks on a 16 3 x 48 lattice. 
The lattice spacing a = 0.1075(7) fm is determined from 
the heavy quark potential with an input ro = 0.49 fm. 
For the strange quark mass, we choose two different val- 
ues m s — 0.080 and 0.100 in the lattice unit. For the 
former, six values of up and down quark masses m u d = 
0.002, 0.015, 0.025, 0.035, 0.050, and 0.080 are taken. For 
the latter, five values m ud = 0.015, 0.025, 0.035, 0.050, 
and 0.100 are used. The smallest value m u d = 0.002 
roughly corresponds to 3 MeV in the physical unit, with 
which pions are in the e regime. In order to investigate 
the finite volume scaling, we also simulate on a 24 3 x 48 
lattice at the same lattice spacing with one choice of the 
sea quark masses m u d = 0.025 and m s — 0.080. 

In the hybrid Monte Carlo (HMC) updates, the global 
topological charge Q of the gauge field is fixed to its initial 
value by introducing extra (unphysical) Wilson fermions, 
which have a mass of cutoff order 16]. In our main runs, 
we set Q = 0. We also simulate another sector of topo- 
logical charge Q = 1 at m u d = 0.015 and m s = 0.080. 

We accumulate 2500 HMC trajectories for the main 
runs in the p regime, 4750 (but the trajectory length is 
0.5) for the e regime lattice, 1800 for the Q = 1 run, 
and 1900 on the 24 3 x 48 lattice. Eighty pairs of low- 
lying eigenvalues of the massless overlap operator D are 
computed at every 5 (or 10 in the Q = 1 and L = 24 
runs) trajectories. For the comparison with ChPT, every 
complex eigenvalue A°" is projected onto the imaginary 
axis as A = lmX ov /(l - Re\ ov / (2m )) • Here m (= 1.6) 
is a parameter to define the overlap-Dirac operator. In 
the analysis, we consider positive A only. The integrated 
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0.002 
0.015 
0.015 
0.025 
0.025 
0.035 
0.050 
0.080 



0.00204(08 
0.00314(18 
0.00354(48 
0.00333(18 
0.00306(07 
0.00404(39 
0.00423(22 
0.00453(23 



0.015 
0.025 
0.035 
0.050 
0.100 



0.00309(14 
0.00349(20 
0.00418(40 
0.00383(13 
0.00520(22 



0.0465(100) 
0.0536(15 
0.0521(25 
0.0624(20 
0.0616(40 
0.0636(17 
0.0696(16 
0.0767(14 



0.0564(19 
0.0622(17 
0.0647(14 
0.0713(16 
0.0835(14 



0.00204(06) 
0.00305(17) 
0.00319(58) 
0.00326(18) 
0.00304(07) 
0.00393(36) 
0.00413(21) 
0.00444(22) 



0.0423(49) 
0.0551(16) 
0.0558(62) 
0.0647(20) 
0.0645(41) 
0.0666(16) 
0.0738(16) 
0.0828(14) 



(Q = i) 

(L=24) 



0.00303(13) 
0.00342(19) 
0.00409(38) 
0.00376(13) 
0.00500(19) 



0.0578(19) 
0.0642(17) 
0.0673(14) 
0.0747(16) 
0.0924(16) 



TABLE I: Numerical results for E e ff and F. The upper half 
is the data at m s = 0.080 while the lower is at 0.100. 



autocorrelation time of the lowest A is measured as 6- 
24 trajectories depending on the simulation parameters. 
The statistical error is estimated by the jackknife method 
after binning data in every 100 trajectories. Details of the 
numerical simulation will be reported elsewhere. 

At each set of sea quark masses, the formula ([TJ is 
described by two unknown quantities S G ff and F. Note 
that X in (3) can be replaced by £ c ff neglecting higher 
order effects. We first determine these parameters from 
the lattice data of pq ( A) . Roughly speaking, the height of 
Pq{X) near A = determines E e g according to the Banks- 
Casher relation, while the shape in the bulk region is 
controlled by F, as far as A is in the region of convergence 
of the chiral expansion. 

Figure [T] shows the spectral density pq(A) multiplied 
by 7r (left panel) and the mode number below A, Nq(X) = 

V J„ A dX'pQ(X') (right), calculated at m U( i — 0.015 and m s 
= 0.080. The solid curve represents the ChPT result ^ 
with S c ff and F determined from Nq(X) at two reference 
points A = 0.004 (~ 7 MeV) and 0.017 (~ 30 MeV). We 
observe that the formula ([TJ) describes the lattice data 
well in the region below A ~ 0.03 (~ m s /2) The result 
is stable within statistical error under changes of the ref- 
erence points in the range A < 0.03. Beyond this value, 
higher order effects may become larger as suggested in 
the analysis of the pion mass and decay constant [13] • In 
the same figure, we also draw the first term of ([T]). Its 
discrepancy from the lattice data for A > 0.01 indicates 
that the second term p p (X, {m sea }) in ([1]) is important 
for the consistency between QCD and ChPT. 

Results from the e regime run are shown in Fig. [2l E e ff 
and F are determined from reference points A — 0.01 and 
0.02. We observe a good agreement between the data and 
NLO ChPT. The NLO correction is less significant than 
that in the p regime, but still visible above A ~ 0.02. 

The values of £ er r and F are summarized in Table U for 
all parameter choices. We use the ChPT formulas of both 
N f = 2 and 2+1 cases. The N f = 2 ChPT formula is 
understood as the leading contribution in an expansion 
in terms of the large strange quark mass. £ and F in 
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FIG. 1: Spectral density npQ(X) (left) and mode number JVg(A) (right) of the Dirac operator at m u d = 0.015, m s — 0.080, and 
Q — 0. The lattice result (given by histogram (left) or solid symbols (right)) is compared with the ChPT formula |T} drawn 
by solid curves. For comparison, the prediction at the leading order of e expansion (dashed curves) is also shown. 
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FIG. 2: Same as the left panel of Fig. [TJ but calculated on 
the e regime lattice, at m u d = 0.002 and m s = 0.080. 



this framework depend on the strange quark mass. The 
curves in Figs. [1] and [2] are drawn using the Nf = 2 + 1 
formula, but the difference from Nf — 2 is hardly visible 
in the range A < 0.03. The numerical results of £ c ff 
and F are, in fact, insensitive to the choice of Nf in the 
formula, except for F in the heavy mass region. We also 
note that there is no significant difference of S c ff between 
m s = 0.080 and 0.100, which confirms decoupling of the 
strange quark from the low-energy dynamics. 

From the data in the non-trivial topological sector 
Q = l,we observe that the topological charge Q largely 
affects the spectral density near A ~ 0, but the values 
of £ e ff and F are consistent with those at Q — 0, as 
listed in Table fl] The data at L = 24 also show the ex- 
pected scaling behavior from ([1]). Since the definition of 
£ c ff O explicitly contains the lattice volume, the results 
from different volumes cannot be compared directly. Af- 
ter converting the L — 24 lattice result E ff = 0.00306(7) 
to that of L — 16, we obtain 0.00341(18), which is con- 
sistent with 0.00333(18) obtained on the L = 16 lattice. 



FIG. 3: Three parameter fit of E eff to the N f = 2 + 1 ChPT. 
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TABLE II: Systematic errors for [E phys (2 GeV)] 1/3 . The total 
error is obtained by adding each estimate by quadrature. 



Next, we analyze the sea quark mass dependence of 
£ e ff from which S, F and Lq can be determined. To see 
the convergence of the chiral expansion, we carry out fits 
using four, five and six lightest data points as a function 
of m u d with m s fixed at 0.080. The data points and fit 
curves of the Nf = 2 + 1 formula are shown in Fig. [3] 
The curvature due to the chiral logarithm in Q is man- 
ifest.The fit result for S? 11 ^, which is S cff in the limit of 
V = oo and m uc i = while keeping m s fixed at 0.08, is 
stable under change of the fitting range. Since we observe 
no sizable m s dependence (see Tabled]), S phys can be con- 
sidered as the one at the physical strange quark mass. 
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From the five points fit, we obtain E phys = 0.00186(10), 
F = 0.0406(5) and L 7 6 = -0.00011(25) in the lattice unit, 
with x 2 / dof = 0-7. 

Since F appears starting at the NLO correction in the 
formula, m u d dependence of the data given in Table U 
reflects the NNLO effects, which is beyond the scope of 
this work. A naive linear extrapolation to the chiral limit 
yields F = 0.0410(46), which roughly agrees with the 
value from the fit of E G ff ■ 

Our final result for the chiral condensate £P h y s , in the 
limit of m u( ]_ — and m s fixed at its physical value, is 

E^(2 GeV) = [242(04)(±i|) MeV] 3 , (5) 

where the errors are statistical and systematic, respec- 
tively. The lattice scale a = 0.1075(7) fm is determined 
from the heavy quark potential ro = 0.49 fm. We use 
the nonperturbatively calculated rcnormalization factor 
l/Z s {2 GeV) = 0.806(12) (ij^) [3 to convert the result 
to the MS scheme at 2 GeV. 

Possible systematic errors are listed in Table HI The 
error from chiral fit is estimated by taking variations of 
the fitting range and the choice of Nf in the ChPT for- 
mula. Finite volume effect is estimated by taking the 
difference between the data on 16 3 x 48 and 24 3 x 48 lat- 
tices. The discretization effect is hard to estimate within 
the calculation done at a single lattice spacing, but partly 
reflected in the mismatch of the lattice spacing obtained 
from different inputs: 0.100(5) fm from the pion decay 
constant [19j and 0.109(2) fm from the f2 baryon mass 
[20j |. To be conservative, the maximum deviation from 
the central value (~ 7.4%) is added in both positive and 



negative directions in Table [TTl 

The chiral condensate obtained in this work ([5]) is 
consistent with other determinations from the pseu- 
doscalar meson mass, £^(2 GeV) = [257(14_)MeV] 3 
and from the top olog ical susceptibility, E MS (2 GeV) = 
[249 (4) MeV] 3 [111, HJ. The former is obtained with the 
NNLO ChPT formula, while the latter only uses the LO 
relation (and the errors do not contain the systematic 
effects). Our result is also consistent with two- flavor re- 
sults in the previous works P, H, [H, [13, HE HI] • Namely, 
there is no significant effect of the strange sea quark. 

We also obtain F = 74(1) (8) MeV and £§(770 MeV) = 
—0.0001 1(25) (11). Their systematic errors are estimated 
in a similar manner. 

By the use of the eigenvalue density of the Dirac op- 
erator calculated on the lattice, the chiral condensate is 
determined without suffering from large subtraction of 
ultraviolet divergences. The dependence on the volume, 
topological charge and quark masses is well described by 
ChPT at NLO in the region where both A and m u d are 
smaller than m s /2. 
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